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A  MODEL  FOR  ASYMMETRIC  HYSTERESIS  IN  PIEZOCERAMIC  MATERIALS 

RALPH  C.  SMITH*  AND  ZOUBEIDA  OUNAIESt 

Abstract.  This  paper  focuses  on  the  characterization  of  hysteresis  exhibited  by  piezoelectric  materials 
at  moderate  to  high  field  levels.  For  soft  materials  in  which  dipoles  are  easily  reconfigured,  the  hysteresis 
loop  is  observed  to  be  rotationally  symmetric  about  the  zero  field,  zero  polarization  point  and  symmetric 
models  can  be  employed.  In  harder  materials,  however,  the  loops  are  no  longer  rotationally  symmetric 
which  necessitates  the  development  of  commensurate  characterization  techniques.  The  model  considered 
here  is  based  upon  the  quantification  of  reversible  and  irreversible  changes  in  polarization  due  to  bending 
and  translation  of  domain  walls  pinned  at  inclusions  inherent  to  the  materials.  The  performance  of  the 
model  is  illustrated  through  comparison  with  PZT4  data. 
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1.  Introduction.  An  inherent  property  of  piezoelectric  materials  is  the  presence  of  hysteresis  and 
constitutive  nonlinearities  at  moderate  to  high  drive  levels.  In  certain  applications,  the  degree  of  hysteresis 
can  be  minimized  by  restricting  input  field  or  stress  levels,  or  employing  certain  feedback  mechanisms  [7]. 
However,  to  attain  the  full  capabilities  of  the  materials  in  high  performance  applications,  it  is  necessary 
to  quantify  the  hysteresis  and  nonlinearities  in  a  manner  which  permits  both  material  characterization  and 
efficient  control  design.  The  latter  criterion  is  facilitated  by  the  development  of  models  which  can  either 
be  exactly  or  approximately  inverted  to  permit  the  construction  of  inverse  compensators  for  linear  control 
design. 

A  common  form  of  hysteresis  in  piezoelectric  materials  is  the  sigmoid  relation  between  the  input  field 
E  and  polarization  P  illustrated  in  Figure  3.1.  For  soft  materials  in  which  dipoles  are  easily  reversed, 
the  curves  are  observed  to  be  approximately  rotationally  symmetric  about  the  point  E  =  0,  P  —  0.  For 
harder  materials  which  are  initially  poled,  the  relation  is  often  asymmetric  due  to  the  fact  that  the  preferred 
orientation  induced  during  poling  is  less  easily  modified  by  AC  fields  at  normal  operating  temperatures. 
Such  a  response  is  illustrated  in  Figure  3.1a  with  data  from  a  poled  PZT4  wafer.  For  general  applications,  it 
is  advantageous  to  have  a  model  which  accommodates  both  the  symmetric  unpoled  response  and  asymmetric 
poled  response  through  the  characterization  of  physical  mechanisms  which  produce  the  hysteresis.  This  is 
accomplished  here  through  the  quantification  of  energy  losses  which  occur  when  domains  reorient  and  domain 
walls  move  in  response  to  an  applied  field.  Details  regarding  the  associated  reversible  and  irreversible  domain 
wall  mechanisms  can  be  found  in  [1,  6]. 

Various  techniques  have  been  employed  for  modeling  hysteresis  in  piezoelectric  materials  including  macro¬ 
scopic  phenomenological  approaches  [3,  4,  15],  microscopic  energy-based  models  at  the  lattice  or  grain  level 
[8],  or  semi-macroscopic  theories  which  employ  energy  relations  in  combination  with  macroscopic  averages 
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to  yield  models  having  effective  parameters  which  characterize  the  bulk  behavior  of  the  material  [2].  The 
present  model  fits  in  this  latter  category. 

The  initial  theory  underlying  this  model  was  developed  in  [10, 11]  for  general  ferroelectric  materials  where 
it  was  illustrated  for  PMN  at  low  temperatures.  This  theory  quantifies  through  electrostatic  energy  relations 
the  reversible  and  irreversible  effects  of  domain  wall  bending  and  translation.  Macroscopic  averages  then 
provide  a  model  having  five  parameters  which  quantifies  the  bulk  attributes  of  the  material.  The  extension 
of  this  model  to  piezoelectric  compounds  and  validation  for  PZT5A  was  provided  in  [12].  The  initial  models 
and  validation  results  accommodate  only  symmetric  hysteresis  loops. 

The  extension  of  the  model  to  accommodate  asymmetric  hysteresis  loops  of  the  type  depicted  in  Figure 
3.1  is  outlined  in  the  next  section  and  then  illustrated  through  a  comparison  with  data  from  both  poled 
and  depoled  PZT4  wafers.  We  note  that  the  basis  of  the  technique  on  physical  mechanisms  yields  a  model 
in  which  the  effects  of  poling  are  accommodated  by  biases  in  the  field  and  polarization  with  the  original 
five  hysteresis  parameters  remaining  unchanged.  This  provides  the  model  with  significant  flexibility  for 
both  material  characterization  and  control  design  since  changes  due  to  poling  (or  depoling)  are  incorporated 
through  the  physics  rather  than  changes  in  the  model  or  hysteresis  parameters. 

2.  Model  Development.  The  hysteresis  between  the  electric  field  E  and  polarization  P  is  modeled  in 
two  steps.  In  the  first,  the  ideal  anhysteretic  polarization,  which  would  result  in  the  absence  of  domain  wall 
pinning,  is  quantified  by  employing  Boltzmann  statistics  to  determine  the  probability  of  dipoles  occupying 
specified  energy  states.  Hysteresis  is  then  incorporated  by  computing  deviations  from  the  anhysteretic  state 
due  the  restriction  of  domain  wall  movement  by  pinning  sites  in  the  material.  This  approach  is  motivated 
by  analogous  theory  for  ferromagnetic  materials  [5].  An  analogous  model  for  symmetric  hysteresis  loops  in 
ferroelectric  materials  is  summarized  in  [10,  11]  with  initial  validation  for  piezoelectric  materials  provided  in 
[12,  13]. 

To  characterize  the  anhysteretic  polarization  which  would  result  in  the  absence  of  pinning  sites,  it  is 
necessary  to  consider  the  effective  field  acting  on  dipoles  in  the  material.  In  the  original  symmetric  models, 
the  effective  field  was  taken  to  be 


Ee  =  E  -f  aP 


where  the  component  aP  quantifies  the  field  contributions  due  to  interdomain  coupling  and  certain  stress 
effects.  As  detailed  in  [10,  11],  the  parameter  a  can  be  formulated  in  terms  of  a  scaling  electric  field  E 
and  the  saturation  polarization  P8 ,  but  is  typically  estimated  in  applications  through  a  least  squares  fit  to 
data.  The  asymmetries  due  to  poling  are  incorporated  by  noting  that  they  produce  a  bias  field  Eq  and  bias 
polarization  Po  which  causes  the  positive  and  negative  coercive  fields  and  remanence  points  to  differ.  The 
resulting  effective  field  is  then  modeled  by 


(2.1) 


Pe  =  (P-P0)  +  a(P-P0). 


As  detailed  in  [10,  11],  the  balancing  of  electrostatic  and  thermal  energy  using  Boltzmann  principles 
yields  a  family  of  models  widely  differ  according  to  assumptions  concerning  possible  domain  orientations.  In 
the  presence  of  a  scaling  polarization  Pi ,  the  assumption  that  cells  can  be  oriented  in  any  direction  yields 
the  Langevin  model 


(2.2) 


Pi+P* 


a 

E~e 
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while  the  assumption  that  dipoles  can  orient  only  in  the  direction  of  the  field,  or  opposite  to  it,  yields  the 
Ising  spin  model 

(2.3)  Pan  =  Pi  +  Pstanh 

In  both  cases,  Ps  denotes  the  saturation  polarization  while  a  is  a  parameter  which  is  estimated  through  a 
least  squares  fit  to  data.  Details  regarding  the  behavior  of  these  two  models  can  be  found  in  [10,  11]. 

Hysteresis  is  incorporated  by  computing  the  irreversible  polarization  P*rr  and  reversible  polarization 
Prev  resulting  from  the  translation  and  bending  of  domain  walls  pinned  at  inclusions  in  the  material.  The 
quantification  of  energy  required  to  break  pinning  sites  yields  the  differential  equation 

/n  a\  dPirr  ~ _ Pan  ~  Pirr _ 

1  j  ~&E  ”  kS  -  a  ( Pan  ~  Pirr) 

specifying  the  irreversible  polarization.  The  parameter  8  =  sign (dE)  ensures  that  the  energy  required  to 
break  pinning  sites  always  opposes  changes  in  polarization  while  5  =  1  if  {dE  >  0  and  P  <  Pan}  or  {dE  < 
0  and  P  >  Pan}  and  5  =  0,  otherwise.  The  parameter  k,  which  quantifies  the  average  energy  required  to 
reorient  domains,  is  demonstrated  in  [13]  to  be  asymptotically  approximated  by  the  coercive  field  Ec  in  soft 
materials. 

The  second  component  of  the  polarization  is  the  reversible  polarization  which  models  the  effects  of 
domain  wall  bending.  To  first  approximation,  this  is  modeled  by  the  relation 

(2.5)  P rev  ~  C(P an  Pirr ) 

where  c  is  a  parameter  which  must  be  estimated  for  the  specific  application. 

The  total  polarization  is  then  given  by 

(2.6)  P  =  Prev  +  Pirr  - 

To  implement  the  model,  the  effective  field  for  a  given  field  and  irreversible  polarization  level  is  computed 
using  (2.1).  This  effective  field  value  is  then  employed  in  either  (2.2)  or  (2.3)  to  compute  the  corresponding 
anhysteretic  polarization.  The  subsequent  irreversible  polarization  is  determined  by  numerically  integrating 

(2.4) .  The  total  polarization  is  then  specified  by  (2.6). 

3.  Model  Validation.  The  polarization  model  (2.6)  provides  the  capability  for  characterizing  hystere¬ 
sis  loops  which  are  either  symmetric  or  asymmetric  functions  of  the  input  field.  To  illustrate  the  performance 
of  the  model,  we  consider  the  characterization  of  both  poled  and  unpoled  PZT4  wafers.  The  reported  results 
were  obtained  from  a  rectangular  wafer  having  a  length  of  3.81  cm,  width  of  0.635  cm  and  thickness  of  0.0381 
cm.  The  wafer  was  initially  poled  in  the  transverse  direction. 

Input  voltages  ranging  from  600  V  to  2200  V  were  applied  to  the  wafer  under  thermally  controlled 
conditions  and  the  resulting  polarization  levels  were  measured.  To  minimize  frequency  effects,  the  input 
frequencies  were  limited  to  200  mHz  and  1  Hz.  While  some  frequency  dependence  was  noted  at  1  Hz,  it 
appears  to  be  manifested  primarily  in  the  field  attained  in  the  material,  and  the  model  performed  adequately 
when  the  field  measured  in  the  material  was  employed  as  input.  We  note  that  under  quasistatic  conditions, 
the  field  generated  in  the  material  by  an  input  voltage  V  can  be  approximated  by  the  relation  E  -V/d  where 
d  =  0.0381  cm  is  the  thickness  of  the  wafer.  The  validity  of  this  relation  degrades,  however,  as  frequency 
effects  are  manifested. 
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Initial  data  sets  were  collected  at  1  Hz  under  the  biased,  and  hence  asymmetric,  conditions  attributed 
to  poling.  The  model  parameters  were  estimated  through  a  least  squares  fit  to  the  2200  V  input  data  which 
yielded  the  values  a  =  6.4  x  106  Vm/C,  a  =  8.0  x  105  C/m2,  k  =  1.5  x  106  C/m2,  c  =  0.5,  Ps  =  -44  C/m2 
and  E0  =  -4.0  x  105  V/m,  P0  =  .01  C/m2,  Pi  =  0.  The  model,  with  these  parameter  values,  was  then  used 
to  predict  the  hysteresis  curves  for  the  material  using  the  measured  field  values  corresponding  to  1000  V, 
1400  V  and  1800  V  inputs.  The  resulting  model  predictions  are  compared  with  the  measured  data  in  Figure 
la.  It  is  observed  that  while  the  model  over-predicts  the  maximum  polarization  for  the  1000  V  input,  it 
characterizes  the  asymmetry  and  accurately  predicts  the  hysteresis  at  the  intermediate  drive  levels. 

The  wafer  was  then  partially  depoled  by  subjecting  it  to  a  sustained  high  voltage  under  thermally 
controlled  conditions  so  that  subsequent  data  was,  at  least  approximately,  rotationally  symmetric  about 
E  =  0,  P  —  0  as  illustrated  in  Figure  lb.  The  model  with  the  same  hysteresis  parameters  a,k,c,  a,Ps 
estimated  from  the  poled  1  Hz,  2200  V  input  data,  was  then  used  with  the  bias  values  E0  =  P0  =  Pi  =  0 
to  predict  the  hysteresis  in  the  depoled  wafer  at  various  input  levels.  The  data  and  model  behavior  for 
the  1000  V  and  1800  V  inputs  can  be  compared  with  those  in  Figure  la,  and  the  lower  voltage  behavior  is 
included  to  further  illustrate  the  flexibility  of  the  model.  We  note  that  while  the  model  fits  in  Figure  lb 
can  be  improved  if  parameters  are  refined,  the  accuracy  attained  using  parameters  estimated  under  different 
operating  conditions  attests  to  the  flexibility  provided  by  the  incorporation  of  associated  physics  in  the 
model. 


4.  Concluding  Remarks.  This  paper  outlines  extensions  to  a  domain  wall  model  for  piezoelectric 
materials  to  accommodate  asymmetries  due  to  poling.  This  modeling  strategy  is  based  on  the  quantification 
of  reversible  and  irreversible  polarization  changes  due  to  domain  wall  pinning  at  inclusions  inherent  to  the 
materials.  The  resulting  model  has  five  hysteresis  parameters  and  three  field  and  polarization  parameters  to 
incorporate  biases  due  to  poling. 

The  performance  of  the  model  was  illustrated  through  a  comparison  with  asymmetric  and  symmetric 
PZT4  data  collected  at  several  drive  levels.  The  model  parameters  were  estimated  through  a  least  squares 
fit  to  data  collected  from  the  poled  wafer  at  a  high  drive  level  and  the  ensuing  model  was  used  to  predict 
both  asymmetric  and  symmetric  hysteresis  at  several  drive  levels.  As  illustrated  by  the  model  fits,  the 
formulation  of  the  model  based  on  energy  principles  provides  it  with  the  capability  for  accurately  predicting 
the  quasistatic  hysteresis  curves  under  varied  poling  and  drive  conditions. 

Finally,  the  ODE  nature  of  the  model  makes  it  amenable  to  inversion  through  the  consideration  of  a 
complementary  ODE  in  a  manner  analogous  to  that  described  in  [9].  This  facilitates  the  construction  of  an 
inverse  compensator  which  can  be  used  for  linear  control  design  [14]. 
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